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Abstract 

My previous calculations of the Lamb shift in muonic hydrogen are reviewed and compared 
with other work. In addition, numerical results for muonic deuterium are presented. 

Introduction 

The energy levels of muonic atoms are very sensitive to effects of quantum electrodynamics 
(QED), nuclear structure, and recoil, since the muon is about 206 times heavier than the 
electron |L. A number of theoretical analyses of the Lamb shift (the 2p-2s transition) 
in light muonic atoms have been published El IH 13 El UH HI M CHI • m °st recently in 
view of a proposed measurement of the Lamb shift im muonic hydrogen J§|. The present 
paper repeats the independent recalculation of some of the most important effects |3] and 
extends the numerical calculations to the case of muonic deuterium, including effects that 
were not considered previously jTUj . Muonic deuterium is in many ways similar to muonic 
hydrogen, but there are some differences. In addition to the different mass the deuteron 
has spin 1 and both magnetic and quadrupole moments. 

In the numerical calculations the fundamental constants from CODATA 2002 
or 1 , he, m M , m e , m u = 137.0359991, 197.32697 MeV fm, 105.658369 MeV, 
0.5109989 MeV, 931.5050 MeV, respectively 

Also, the following properties of the proton and deuteron were used: m p = 938.272 MeV/ c 2 , 
R p = 0.875 ± 0.007 fm (other recent values are discussed below) and fi p = 2.79285 hn- 
Also, m d = 1875.613 MeV/c 2 , R d = 2.139 ± 0.003 fm and fx d = 0.85744 p N = 0.307012 /x p . 
([HI) The deuteron has spin 1 and thus has both magnetic and quadrupole moments. 
The quadrupole moment of the deuteron is taken to be Q = 0.2860(15) fm 2 [T21 HSl El- 



Vacuum Polarization 

The most important QED effect for muonic atoms is the virtual production and anni- 
hilation of a single e + e~ pair It has as a consequence an effective interaction of order 
aZa which is usually called the Uehling potential (|T71 This interaction describes 
the most important modification of Coulomb's law. Numerically it is so important that it 
should not be treated using perturbation theory; instead the Uehling potential should be 
added to the nuclear electrostatic potential before solving the Dirac equation. However, 
a perturbative treatment is also useful in the case of very light atoms, such as hydrogen. 



However, unlike some other authors, we prefer to use relativistic (Dirac) wave functions 
to describe the muonic orbit. Since these contributions have been extensively discussed 
in the literature PjJ |2J El H] (among others), there is no need to go into detail here. The 
results, calculated as the expectation value of the Uehling potential using point-Coulomb 
Dirac wave functions with reduced mass are, for muonic deuterium: 





point nucleus 


R d =2.139 fm 


Uehling 
Kaellen-Sabry 


2px/ 2 - 2s 1/2 2p 3/2 ~ 2Si/ 2 
227.6577 227.6635 
1.66622 1.66626 


2p 1/2 - 2sy 2 2p 3/2 - 2si /2 
227.5985 227.6043 
1.66577 1.66582 



The effect of finite proton size calculated here can be parametrized as -0.0129(r 2 ). However 
higher iterations can change these results. Up to now, these have not been calculated well 
for muonic deterium, as far as I know. 

Corresponding numbers for muonic hydrogen, calculated as the expectation value of 
the Uehling potential using point-Coulomb Dirac wave functions with reduced mass are: 





point nucleus 


J R p =0.875fm 


Uehling 
Kaellen-Sabry 


2pi /2 - 2si /2 2p 3/2 - 2si /2 
205.0282 205.0332 
1.50814 1.50818 


2p 1/2 - 2sy 2 2p 3/2 - 2s 1/2 
205.0199 205.0250 
1.50807 1.50811 



The effect of finite proton size calculated here can be parametrized as -0.0109(r 2 ). How- 
ever higher iterations can change these results. The contribution due to two and three 
iterations have been calculated by 4J and |23j . respectively, giving a total of 0.151 meV. 
An additional higher iteration including finite size and vacuum polarization is given in 
ref. jl] (equations(66) and (67)) and ref. [2| (equations(264) and (268)). These amount to 
-0.0164(r 2 ). The best way to calculate this would be an accurate numerical solution of 
the Dirac equation in the combined Coulomb-plus Uehling potential. 

The mixed muon-electron vacuum polarization correction ([SUE]) is 0.00007 meV for 
hydrogen and 0.00008 meV for deuterium. 

The Wichmann-Kroll contribution was calculated using the parametrization for the 
potential given in pQ. The result obtained for hydrogen is -0.00103 meV, consistent with 
that given in [2|. For deuterium, the contribution is -0.00111 meV. 

The equivalent potential for the virtual Delbriick effect was recomputed from the 
Fourier transform given in [22] arid PQ. The resulting potential was checked by reproducing 
previously calculated results for the 2s-2p transition in muonic Helium, and the 3d-2p 
transitions in muonic Mg and Si. The result for hydrogen is +(0.00135 ± 0.00015) meV, 
and for deuterium it is +(0.00147 ± 0.00016) meV. As in the case of muonic helium, 
this contribution very nearly cancels the Wichmann-Kroll contribution. The contribution 
corresponding to three photons to the muon and one to the proton should be analogous 
to the light by light contribution to the muon anomalous moment; to my knowledge, 
the corresponding contribution to the muon form factor has never been calculated. It 
will be comparable to the other light by light contributions. This graph was included in 
contributions to the muon's anomalous magnetic moment; the contribution to the muon 
form factor is one of the most significant unknown corrections. 
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The sixth order vacuum polarization corrections to the Lamb shift in muonic hydrogen 
have been calculated by Kinoshita and Nio "2H|. Their result for the 2p-2s transition (in 
hydrogen) is 

AE (G) = 0.120045- (aZf-m r « 0.00761 meV 

and 0.00804 meV for muonic deuterium. 

However, I should remark that the contributions from figures 1 and 2 of Ref. j2S] were 
checked by direct integration. Although the results agreed perfectly for the case of hydro- 
gen, there were small, but significant discrepancies for the case of deuterium, (hydrogen: 
Fig. 1 contributes 0.000396 meV and Fig. 2 contributes 0.002931 meV; deuterium: direct 
integration gave 0.000472 meV and 0.003364 meV, respectively, while the work of ref. [23] 
indicates values 0.000419 meV and 0.003906 meV, respectively). This indicates that, at 
least for these two graphs, integration over momentum transfer involves more than a single 
reduced mass factor. 

The hadronic vacuum polarization contribution has been estimated by a number 
of authors j23 I2H 12]- It amounts to about 0.012 meV in hydrogen and 0.013 meV in 
deuterium. One point that should not be forgotten about the hadronic VP correction is 
the fact that the sum rule or dispersion relation that everyone (including myself) used does 
not take into account the fact that the proton (nucleus) can in principle interact strongly 
with the hadrons in the virtual hadron loop. It is irrelevant for the anomalous magnetic 
moment but probably not for muonic atoms. An estimation of this effect appears to be 
extremely difficult, and could easily change the correction by up to 50%. Eides et al. [2] 
point out that the graph related to hadronic vacuum polarization can also contriibute 
to the measured value of the nuclear charge distribution (and polarizability). It is not 
easy to determine where the contribution should be assigned. This may also be true for 
the so-called " proton self-energy" jSJ [2] , which involves some of the same graphs as are 
present in the calculation of radiative corrections to electron scattering. 



Finite nuclear size and nuclear polarization 

The main contribution due to finite nuclear size has been given analytically to order (aZ) 6 
by Friar [21]. The main result is 



AE ns 



2aZ ( aZm r 



n 



(r 2 ) - ^(r 3 ) (2) + (aZ)\E REL + m 2 r F NR ) 



where (r 2 ) is the mean square radius of the proton. For muonic hydrogen, the coef- 
ficient of (r 2 ) is 5.1975 (meV fm -2 ), giving an energy shift (for the leading term) of 
(3.979±0.076) meV if the proton rms radius is 0.875 fm. Other values of the proton radius 
that have been reported recently in the literature are 0.880 fm [2H] and (0.895 ± 0.018 fm) 
[2E]. The second term in Eq.Q contributes -0.0232 meV for a dipole form factor and 
-0.0212 meV for a Gaussian form factor. The parameters were fitted to the proton rms 
radius. This can be written as -0.0347(r 2 ) 3 / 2 or -0.0317(r 2 ) 3 / 2 . This differs slightly from 
the value given by Pachucki [5 j. The model dependence introduces an uncertainty about 
±0.002 meV. The remaining terms contribute 0.00046 meV. This estimate includes all of 
the terms given in [2*4] . while other authors [5J give only some of them. Clearly the 
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neglected terms are not negligible. There is also a contribution of -3 • 10~ 6 meV to the 
binding energy of the 2p 1 / 2 -level, and a recoil correction of 0.013 meV to the binding 
energy of the 2s-level. 

Pachucki j3] has estimated a correction similar to the second term (proportional to 
(r 3 )( 2 )) in Eq.(fTj). Since the logarithmic terms in the two-photon correction without finite 
size (see below) also seem to be suspect, this correction requires further investigation. 
In particular, the parametrization of the form factors used in any calculation should 
reproduce the correct proton radius. 

For muonic deuterium, the main contribution amounts to 
-6.0732 (r 2 ) =-(27. 787±0. 078) meV. Depending on the model, the term proportional to 
(r 3 )(2) gives a contribution of 0.382 meV or 0.417meV. 

As mentioned previously, the finite-size contributions to vacuum polarization in muonic 
hydrogen can be parametrized as — 0.0109(r 2 ) — 0.0164(r 2 ), giving a total of — 0.0273(r 2 ) 
or -0.0209(6) meV if the proton radius is 0.875 fm. For deuterium, only the contribution 
corresponding to the first term of the sum (— 0.0129(r 2 )) has been calculated. 

The contribution due to nuclear polarization (in hydrogen) has been calculated by 
Rosenfelder [H to be 0.017 ± 0.004 meV, and by Pachuki to be 0.012 ± 0.002 meV. 
Other calculations [SHI Ell §i ye intermediate values (0.013 meV and 0.016 meV, respec- 
tively). The value appearing in tabled is an average of the three most recent values, with 
the largest quoted uncertainty, which is probably underestimated. 



Relativistic Recoil 

As is well-known, the center-of-mass motion can be separated exactly from the relative 
motion only in the nonrelativistic limit. Relativistic corrections have been studied by 
many authors, and will not be reviewed here. The relativistic recoil corrections summa- 
rized in P include the effect of finite nuclear size to leading order in m^/rriN properly. 

Up to now this method has been used to treat recoil corrections to vacuum polarization 
only in the context of extensive numerical calculations that include the Uehling potential 
in the complete potential, as described in [I]. They can be included explicitly, as a 
perturbation correction to point-Coulomb values. Recall that (to leading order in 1/mjy), 
the energy levels are given by 

£ = £ '-2^ + 2^<' ! ( r » + 2B » P ' M > < 2) 

where E r is the energy level calculated using the reduced mass and B is the unperturbed 
binding energy. Also 

h(r) = -Pi(r)(Pi(r) + -Q 2 (r)) - Ig^r^r) + Q 4 (r)/r 3 ] (3) 
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Here 



P^r) = 4naZ r'p(r')dr' = - V (r) - rV '(r) (4) 

J T 

rr 

g 2 (r) = AnaZ / r' 2 p{r')dr' = rV(r) 

Jo 

rr 

Q 4 (r) = AnaZ / r' A p{r')dr 
Jo 

An effective charge density p^p f° r vacuum polarization can be derived from the 
Fourier transform of the Uehling potential. Recall that (for a point nucleus) 

V VM (r) = -~-X 1 (2m.r) 



3tt J x z 2 \ 2z 2 J \n J q 2 + Am 2 z 2 

where Xn{%) is defined in [I]. In momentum space, the Fourier transform of V 2 V is 
obtained by multiplying the Fourier transform of V by — q 2 . Note that using the nor- 
malizations of [HE!) one nas V 2 V^ = — AiraZp where p is the charge density. One then 
obtains 



2a r j C* 2 - !) 1/2 L 1 W 2 r g 4 -Jo(^) , 



Z ' ,2 



7T 



(5) 



<7 U 2 (q)jo(qr) dq 



[/2(g) is defined in [I]. 

Keeping only the Coulomb and Uehling potentials, one finds 

2a 

Pi(r) = (2m e )x (2m e r) 

37T 

/ 2a 

Q 2 (r) = aZ ( 1 + — [xi(2m e r) + (2m e r)xo(2m e r)] 



2a r , (z 2 - I) 1 ' 2 ( 1 
2\ f°° 1 (6gr — (gr) 3 ) cos(gr) + (3(qr) 2 — 6) sin(gr) 



71 J Jo 9 + ^m e z <? 



dq 



Details of the calculations for the case of vacuum polarization are given in Appendix 1 
and in Ref.jSj. Corrections due to finite nuclear size can be included when a model for 
the charge distribution is given. This done by Friar [21] (and confirmed independently 
for two different model charge distributions); the contribution due to finite nuclear size to 
the recoil correction for the binding energy of the 2s-level is -0.013 meV. The factor l/m n 
is replaced by l/(m^ + m N ), also consistent with the calculations presented in [23] . 

Combining the expectation values given in Appendix 1 according to equations Eland El 
one finds a contribution to the 2p-2s transition of -0.00419 meV (hydrogen) and -0.00479 meV 
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(deuterium). To obtain the full relativistic and recoil corrections, one must add the 
difference between the expectation values of the Uehling potential calculated with rel- 
ativistic and nonrelativistic wave functions, giving a total correction of 0.0166 meV for 
muonic hydrogen. This is in quite good agreement with the correction of .0169 meV 
calculated by Veitia and Pachucki . The treatment presented here has the advantage of 
avoiding second order perturbation theory. For deuterium, one obtains a total correction 
of 0.0179 meV. 

The review by Eides et al. [2] gives a better version of the two photon recoil (Eq. 136) 
than was available for the review by Borie and G. Rinker [T|. Evaluating this expression 
for muonic hydrogen gives a contribution of -0.04497 meV to the 2p-2s transition in 
hydrogen and -0.02656 meV in deuterium. Higher order radiative recoil corrections give 
an additional contribution (in hydrogen) of -0.0096 meV [2j- However, some of the 
contributions to the expressions given in [2] involve logarithms of the mass ratio m^/rriN- 
Logarithms can only arise in integrations in the region from to m^; in this region 
the effect of the nuclear form factor should not be neglected. Pachucki ^ has estimated 
a finite size correction to this of about 0.02 meV, which seems to be similar to the term 
proportional to (r 3 )( 2 ) given in Eq.Q as calculated in the external field approximation 
by Friar [21]. This two-photon correction requires further investigation. In particular, 
the parametrization of the form factors used in any calculation should reproduce the 
correct proton radius. Also the relationship among the different contributions needs to 
be specified more clearly. 



An additional recoil correction for states with t ^ has been given by [33] (see also 
|21). It is 



When evaluated for the 2p-states of muonic hydrogen, one finds a contribution to the 
2p-2s transition energy of 0.0575 meV for the 2piy 2 state and -0.0287 meV for the 2p 3 / 2 
state in hydrogen ( 0.0168 meV for the 2pi/ 2 state and -0.0084 meV for the 2p 3 / 2 state in 
deuterium) 

A final point about recoil corrections is that in the case of light muonic atoms, the mass 
ratio m^/m^v is considerably larger than the usual perturbation expansion parameter aZ. 
Contributions of higher order in the mass ratio could be significant. 

Muon Lamb Shift 

For the calculation of muon self-energy and vacuum polarization, the lowest order (one- 
loop approximation) contribution is well-known, at least in perturbation theory. Including 
also muon vacuum polarization (0.0168 me V) and an extra term of order (Za) 5 as given 
in j2]: which contributes -0.00443 meV, one finds a contribution of -0.66788 meV for the 
2si/ 2 — 2pi/2 transition and -0.65031 meV for the 2si/ 2 — 2p 3 / 2 transition. For deuterium, 
the corresponding contributions are given by -0.77462 meV for the 2sx/2 — %Pi/2 transi- 
tion and -0.75512 meV for the 2sx/2 — 2p 3 / 2 transition. The second order calculation in 
deuterium includes muonic vacuum polarization (0.01968 meV); the extra term of order 
(Za) 5 as given in j2], contributes -0.00518 meV. 
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These results, and the higher order corrections PjJ |^] can be summarized as 



Transition 


2px /2 - 2Sx/2 


2^3/2 - 2si/ 2 


Hydrogen 
second order 
higher orders 


-0.66788 
-0.00172 


-0.65031 
-0.00165 


Total 


-0.66960 


-0.65196 


Deuterium 
second order 
higher orders 


-0.774616 
-0.002001 


-0.755125 
-0.001926 


Total 


-0.776617 


-0.757051 



Table 1: Contributions to the muon Lamb shift (i?(2p 1/ / 2 ) — £J(2s 1( / 2 )) in muonic 
hydrogen and deuterium, in meV. 



For hydrogen, Pachuki [I] has estimated an additional contribution of -0.005 meV for 
a contribution corresponding to a vacuum polarization insert in the external photon. 
The higher order contributions can be written in the form 

AEls = 1_ . (vv > [<m + 1] + ° fm l(l d li • a,) 

where ^(0) = a^; the higher order contributions (fourth and sixth) can be taken from 
the well-known theory of the muon's anomalous magnetic moment: 
F 2 (0) = a ll = a/2n + 0.7658(a/7r) 2 + 24.05(a/7r) 3 . 
The fourth order contribution to ^(0) is 

0.46994(a/vr) 2 + 2.21656(a/7r) 2 = 2.68650(a/vr) 2 [T]. The sixth order contributions to 
^{(0) that involve electron vacuum polarization loops (especially the light-by-light graph) 
might contribute at an experimentally significant level, but have not been calculated. 
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Summary of contributions for muonic hydrogen 

Using the fundamental constants from the CODATA 2002 (jllj) one finds the transition 
energies in meV in table [21 Here the main vacuum polarization contributions are given 
for a point nucleus, using the Dirac equation with reduced mass. Some uncertainties have 
been increased from the values given by the authors, as discussed in the text. 

The finite size corrections for hydrogen up to order (aZ) 5 can be parametrized as 
5.1975(r 2 ) + 0.0109(r 2 ) + 0.0164(r 2 ) + 0.0347(r 3 ) (2) . The various contributions are dis- 
cussed in the text. 



(!nntTirint~inn 


\f. t\ 1 n p ( m p V 1 

V CXI LLC I lllv; V J 


T Tn pprt a intv ( ttipV i 
uiik/Ci ucxiiiuy i iiic v i 


1 oh in rf 
U cllllllg, 


90^ 0989 




xVclllcIl- O a, L)I y 


l.OUo 1 




VV lCIlIIldllll-lVI Oil 


-U.UU1UO 




AAirf TjplnTnppK 
v ii \j . uciui Lie v. i\- 


001 35 

W . WW ± *J W 


0001 5 

W . WWW ± <J 


IIllXt!G IllU-t: V r 


n nnnn7 

U.UUUU ( 




npiHrnnip VP 


01 1 

W . W -L -L 


002 


sixth order [23J 


0.00761 




recoil [2] (eql36) 


-0.04497 




recoil, higher order [2] 


-0.0096 




recoil finite size 24 


0.013 


0.001 


recoil correction to VP jT] 


-0.0042 




additional recoil EI] 


0.0575 




muon Lamb shift 






second order 


-0.66788 




fourth order 


-0.00169 




nuclear size (i? p =0.875fm) 




0.007fm 


main correction [2*1 


-3.979 


0.076 


order (aZf [21 


0.0232 


0.002 


order {aZf [21 


-0.0005 




correction to VP 


-0.0083 




polarization 


0.015 


0.004 


Other (not checked) 

VP iterations H 

VP insertion in self energy H 

additional size for VP [2] 


0.151 
-0.005 
-0.0128 





Table 2: Contributions to the muonic hydrogen Lamb shift (the 2s 1/2 — 2p±/2 
transition). The proton radius is taken from 



Summary of contributions for muonic deuterium 

For deuterium, one finds the transition energies in meV in table [3] Also here the main 
vacuum polarization contributions are given for a point nucleus, using the Dirac equation 
with reduced mass. The finite size corrections for deuterium up to order (aZ) 5 can be 
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parametrized as 

6.0732(r 2 ) + 0.0129(r 2 ) + 0.0409 (r 3 ) (2) , although not all contributions to the effect of 
finite size on the vacuum polarization correction are included. 



Contribution 


Value (meV) 


Uncertainty (meV) 


T T 1 1 * 

Uehlmg 


227.6577 




Kallen-Sabry 


l.OOOZ 




Wichmann-Kroll 


n nm 1 1 
-0.00111 




virt. Delbrueck 


0.00147 


n nnni ^ 

0.00016 


mixeu mu-e v r 


n nnnns 

U.UUUUo 




hadronic VP 


U.Ulo 


U.UUz 


sixth order [23J 


U.UUoU4 




recoil |2J (eql36j 


-0.02656 




recoil, higher order jz] 


? 




recoil, finite size |24j 


0.019 


n nno 
0.003 


lecoii conecLion iu vr jij 


-u.uuio 




annii'inTia Ttir'/^n 1 < /I 
cLLlLllLlUllcll ICUUll IQ^tl 


U.U1UO 




muon Lamb shift 






second order 


-0. 1 /462 




luuiLii oiuei 


n nn9nn 




nuclear size (_R,i=2.139 fm) 




0.003 fm 


main correction [21 


-27.787 


0.078 


order (aZf [21 


0.0400 


0.018 


order {aZf [21 


-0.0045 




correction to VP 


-0.0592 




polarization 


? 




Other (not checked) 






VP iterations d 


? 




VP insertion in self energy jl 


? 




additional size for VP [2] 


? 





Table 3: Contributions to the muonic deuterium Lamb shift. Lhe deuteron radius is 
taken from |llj . 



Fine structure of the 2p state 

The fine structure of the 2p states can be calculated by using the relativistic Dirac 
energies, computing the vacuum polarization contributions with Dirac wave functions, 
and including the effect of the anomalous magnetic moment in the muon Lamb shift. An 
additional recoil correction (EqEJ) also has to be included. The results are summarized 
in table [U One should also include the 5 2 /2Mjv-type correction to the fine structure, 
(see j2j, Eq(38)). This is tiny (5.7- 10~ 6 meV in hydrogen) and is not included in the table. 
Friar [21 has given expressions for the energy shifts of the 2p-states due to finite nuclear 
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Hydrogen 


Deuterium 


Dirac 


o. 41004 


o.o04oU 


uemmg^ v r j 


u.uuou 


U.UUO 1 


Kallen-Sabry 


0.00004 


0.00005 


anomalous moment a M 






second order 


0.01757 


0.01491 


higher orders 


0.00007 


0.00007 


Recoil (Eq.®) 


-0.0862 


-0.0252 


Total Fine Structure 


8.352 


8.864 



Table 4: Contributions to the fine structure (E(2p^/ 2 ) — E(2pi/ 2 )) of the 2p-state in 
muonic hydrogen and deuterium, in meV. 



size. These were calculated and found to give a negligible contribution (3.1 • 10 6 meV) 
to the fine structure of the 2p-state in hydrogen. 



Hyperfine structure 

The Breit equation |3*4*1 12 U\ contributions to the fine- and hyperfine interactions for 
general potentials and arbitrary spins were given by Metzner and Pilkuhn [36J. Here a 



version applicable to the case of muonic atoms {Z\ = —1, s\ 
Z 2 = Z) is given. 



1/2, mi = m„, «i 



1 1 dV 



L.si 



2m fl r dr 



1 + a, 



1 

m,, 



L ■ s i 



(7) 



This can be rearranged to give the well-known form for spin 1/2 particles with an 
anomalous magnetic moment, namely 



1 dV 1 + a M + (a M + l/2)m N /m fM 



r dr 



Note that 



+ 



1 1 



L ■ a„ 



1 



2m?, 



2ml 



2m% 



m N m^ 

so that the terms not involving a M in the spin-orbit contribution are really the Dirac fine 
structure plus the Barker-Glover correction (Eq. EJ). 
Also 



Vi 



1 1 dV 



L,S2 



2m 2 r dr 



1 + K 2 /Z 



s 2 m r 



1 

m 2 



L ■ s 2 



Usually one writes 



Z + K 2 



m r . 



m 2 ,,u P 

where \i 2 is the magnetic moment of the nucleus in units of nuclear magnetons (///v = e/2m p ) 
A value of fid = 0.85744 = 0.307012 /x p corresponds to Kd = 0.714. 



V f 



2(1 + oJ^ 2 



«1,S2 



2s 2 m l] m 2 Yr dr 



1 oV 



(3si • rs 2 • r - si • s 2 ) - ^ VVs\ • s 2 
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r dr 



with Q in units of 1/m 2 ,. The quadrupole moment of the deuteron is taken to be 
Q = 0.2860(15) fm 2 OUnillll • In other units, one finds Q = 25.84/m 2 = 7.345 x 10~ 6 MeV" 2 . 



Note that Vl, s1 describes the fine structure, while the hyperfine structure is described 
(in perturbation theory) by the expectation values of Vl, S2 i V S1>S2 , and Vq (where appli- 
cable) . 

The Uehling potential has to be included in the potential V(r). For states with £ > 
in light atoms, and neglecting the effect of finite nuclear size, we may take 



ldV 

r dr 



aZ 



2a r(z 2 -lW 2 



1 + T" I 

Sir J 1 



1 + p) • (l + 2m e rz)-e 



-2m e rz 



dz 



which is obtained from the Uehling potential fTJ ^] by differentiation. Then, assuming 
that it is sufficient to use nonrelativistic point Coulomb wave functions for the 2p state, 
one finds 

1 \ / 1 



r 3 / 2 P 



r 6 / 2 P 



where 



2a 

e 2p — ~ 



oo ^2 _ ^1/2 



3tt J l z 2 \ 1 + 2z 2 J \(l + az) 2 ' (1 + azf 

with a = 2m e /(aZm r ). For hydrogen, e 2p = 0.000365, and for deuterium e 2p = 0.000391. 



1 



2az 



dz 



(9) 



Hyperfine structure of the 2p state in muonic hydrogen 

The hyperfine structure of muonic hydrogen is calculated in the same way as was done in 
earlier work [3E3, but with improved accuracy. Most of the formalism and results are 
similar to those given by |lj and 

The hyperfine structure of the 2p-state is given by [3 ES] (F is the total angular 
momentum of the state) 



ldV\ 
Am^m 7v \ r dr 1 2p 



■1 + k) 



2(1 + Xp ) Sjj ,(F(F + l)- 11/4) 



+ 6jf(C Fl (l + a„)-2(l + x))l. ! ! 



F 1 



2 2 



t F 1 

2 2 J 



(10) 



where j = y/2j + 1, the 6-j symbols are defined in [SZj , and 
Cpi = 6 fi - 25 FQ - {l/h)8 F 2 
Also 

m M (l + 2k p ) 



0.09245 



2m p (l + k p ) 

represents a recoil correction due to Thomas precession EH ES] • The correction due to 
vacuum polarization (Eq. (JUJ)) should be applied to the HFS shifts of the 2p-states. 
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As has been known for a long time [7J [H3 El ESJ, the states with total angular 
momentum F = 1 are a superposition of the states with j = 1/2 and j = 3/2. Let 
the fine structure splitting be denoted by 5 = E 2p3 / 2 — E 2p i/2 = 8.352 meV, and let 

(a£)V 1 + 

and j3' = {3 P ■ (1 + £2p)- The matrix elements for the hyperfine structure of the 2p-state 
are then given by 



j j' Energy 

1/2 1/2 (/5'/8)(2 + a; p + ^)[-5 Fi o + l/3 5 F , 1 ] 

3/2 3/2 5 + 03'/4) (4 + hx d -a„)[- 1/12^ + 1/20^] 

3/2 1/2 ((3'/24)(l + 2x p -a ti )[V2 5 F , l } 

Then for the 2p- level with j = j' = 1/2 and F — 0, the energy shift is given by 
- (P'/8)(2 + x p + a M ) =-5.971 meV, and for the 2p-level with j = / = 3/2 and F = 2, 
the energy shift is given by 5 + (/3'/80)(4 + 5x p — a M ) = 9.6243 meV. 

For the 2p-levels with F = 1 the corresponding matrix has to be diagonalized. The 
resulting numerical values for the eigenvalues are (A ± R)/2 = 1.846 meV and 6.376 meV, 
where 

A = 5 - P'(x p - a„)/16 
R 2 = [5- (3'(1 + 7x p /8 + a fl /8)/6} 2 + ((3') 2 (1 + 2x p - a M ) 2 /288 



Hyperfine structure of the 2p-state in muonic deuterium 

For the 2p state, the matrix elements of the magnetic hyperfine structure have been given 
by Brodsky and Parsons [S3]. For hydrogen they are the same as those calculated above. 
Here the Uehling potential will be included in the expectation value of 

1_ dV\ 

r dr I 

as discussed above. 
Let 

= 16(1 + /td) a /(f+l)(2£+l) = ^(«Zm r ) 3 = 4.0906meV^ 

(for a point Coulomb potential) 

The matrix elements for the magnetic hyperfine structure are then given by 



j j' Energy 

1/2 1/2 {P D /6){2 + x d + aJ[-6 Ftl/2 + l/2 6 F<3/2 \ 

3/2 3/2 8 + (/3 D /4)(4 + 5x d -a A1 )[-l/65 Fjl/2 -l/15^3/2 + l/105 F) 5/ 2 ] 

3/2 1/2 ((3 D /48)(l + 2x d -a tM )[V25 F , l/2 -V58 F , 3/2 } 
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where x d = (mj;/m d m r ) («<z/(l + K d )) = 0.0248. 

For the evaluation of the contributions of the quadrupole HFS, let 



n/ 1 dV \ 

For a point Coulomb potential, and the 2p-state, €q = aQ(Zam r ) 3 /24 = 0.43243 meV. 
The quadrupole interaction results in energy shifts of 



j f Energy 

1/2 1/2 

3/2 3/2 e Q [5 F>1/2 -4/5 5 Fj3/2 +l/5 5 F>5/2 \ 

3/2 1/2 eq[y/2^ 1/2 -l/v/5^3/2] 

As mentioned before, the Uehling potential has to be included in the potential V(r). 
For states with £ > in light atoms, this can be taken into account by multiplying (3d 
and €q by (1+ e 2p ) where e 2p is given by Eq.©. With a numerical value of e 2p = 0.000391 
for muonic deuterium, the value of e<g is increased to 0.43440 meV and the value of (3d is 
increased to /?^=4.0922 meV. 

Then for the 2p- level with j = j' = 3/2 and F = 5/2, the energy shift is given by 
5 + e Q /5 + ifl' D /m) (4 + 5x d - a„) = 9.373 meV. 

For the 2p-levels with F — 1/2 and F = 3/2, the corresponding matrices have to be 
diagonalized. The resulting numerical values for the eigenvalues are, 
for F = 1/2, -1.3834 meV and 8.5974 meV; for F = 3/2 they are 0.6856 meV and 
8.2410 meV. 



Hyperfine structure of the 2s-state: 

The expectation value of V SlS2 in an ns state with j = 1/2 is 



2a2a(aZ) 3 ml 



When s 2 = 1/2, and [i 2 jm v = (1 + K2)/m2, this reproduces the well-known result for 
muonic hydrogen: 

AE ns = 8 j aZ rm 3 r (1 + ^ (1 + } = (8/n 3 )/3 . (1 + } = (8/n 3 } x 22.8332 meV^ 

(see, for example [2], Eq. (271,277)). The numerical value was calculated for hydrogen. 
For deuterium, with s 2 = 1, the corresponding hyperfine splitting is 

A = 2 }°f^ mr -(l+K d )-(l+Oa)-[F(F+l)-ll/4] = (8/n 3 )x2.04766meKx[F(F+l)-l 
6n- i m fl m 2 , 

for a total splitting of 6.14298 meV in muonic deuterium. This is in reasonably good 
agreement with the result given by Carboni [TU] . 
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As was shown in |2|, the energy shift of the 2s-state in muonic hydrogen is given 

by: 

AE 2s = (3 ■ (1 + a M ) • (1 + e VP + e vertex + e Breit + e F s,rec) • [hi - 38 F0 ]/4 (11) 

The corrections due to QED effects, nuclear size and recoil are analogous for muonic 
deuterium. 

The QED corrections have been discussed by Borie [31 El UH] (see also |3Sj), and are 
given in Appendix 2. 

The correction due to finite size and recoil have been given in |I| as -0.145 meV, while 
a value of -0.152 meV is given in [32]- Ref. [3] also gives a correction as calculated by 
Zemach (JH3j) equal to -0.183 me V, but claims that this correction does not treat recoil 
properly. The Zemach correction is equal to 

e Z em = -2aZm r (r) {2 ) 

where (r)(2) is given in 01231 HI]- Using the value (r)™ = 1.086 ± 0.012 fm from |4*T| . 
gives Ezem = —0.00702, and a contribution of of -0.1742 meV to the hyperfine splitting 
of the 2s state. Including this, but not other recoil corrections to the hyperfine structure 
of the 2s-state gives a total splitting of 22.7806 meV. Additional higher order corrections 
calculated in Ref. (321 amount to a total of -0.0003 meV and are not included here. 

It would be very desirable to understand the reasons for the discrepancy between 
references 0] and |32] i n the calculations of this effect. Also, since the Zemach radius 
seems to be sensitive to details of the electric and magnetic charge distributions [31], 
evaluations performed with a dipole-type form factor may not be good enough. This 
point requires further invesigation. 

For muonic deuterium, the coefficient of (r)(2) is -0.007398 fm -1 , giving, with 
(r) (2 ) = 2.593 ± 0.016 fm from gj, e Zem = -0.01918 ± 0.00012. 

The total hyperfine splitting of the 2s-state of muonic deuterium, including all correc- 
tions, is 

3 

AE 2s = -(3 D ■ (1 + M ) • (1 + S VP + e vert ex + ^Breit + £FS,rec) = 6.0582 meV 
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Transition 


Energy shift in meV 


1 pi/2 - 1 S1/2 


11 11/1 
11.114 


3 Pl/2 - 1 S1/2 


18.931 


3 P3/2 - 1 S1/2 


23.461 


Vl/2 - 3 S1/2 


-11.666 


3 Pl/2 - 3 S1/2 


-3.849 


3 P3/2 - 3 S1/2 


0.681 


5 P3/2 - 3 S1/2 


3.929 



Table 5: Fine- and hyperfine contributions to the Lamb shift in muonic hydrogen. 



Transition 


Energy shift in meV 


2 Pl/2 ~ 2 Sl/2 


2.655 


2 P3/2 ~ 2 S 1/2 


12.636 


Vl/2 - 2 Si/2 


4.724 


4 P3/2 - 2 Sl/2 


12.280 


2 Pl/2 - 4 Sl/2 


-3.403 


2 P3/2 - 4 Sl/2 


6.578 


V/2 - 4 Sl/2 


-1.334 


6 4 
P3/2 - Si/a 


6.222 


6 4 
P3/2 - Si/2 


7.354 



Table 6: Fine- and hyperfine contributions to the Lamb shift in muonic deuterium. 

TablesElandElgive the contributions to the transition energies due to fine and hyperfine 
structure. 

Summary of contributions and Conclusions 

The most important contributions to the Lamb shift in muonic hydrogen, including 
hyperfine structure, have been independently recalculated. A new calculation of some 
terms that were omitted in the most recent literature, such as the virtual Delbruck effect 
|2*2~] and an alternative calculation of the relativistic recoil correction have been presented. 

Numerically the results given in table El add up to a total correction of 
(206.032(6) - 5.225 (r 2 ) + 0.0347 (r 2 ) 3 / 2 ) meV = 202.055±0. 12 meV. (for the value of the 
proton radius from JI]). As is well known, most of the uncertainty arises from the 
uncertainty in the proton radius. 

Numerical results were also given for muonic deuterium. The total correction is 
(228.573(6) - 6.086 (r 2 ) + 0.0409 (r 2 ) 3 / 2 ) meV = 200. 767±0.09meV. The complete depen- 
dence on the deuteron radius is uncertain since contributions from iteration of the po- 
tential are not included. Also, some other contributions are not included, as indicated in 
table H 
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Appendix 1: Details of the Relativistic Recoil Calculation 

As mentioned above, the energy levels of muonic atoms are given, to leading order in 1/mjv by 

where E r is the energy level calculated using the reduced mass and Bq is the unperturbed 
binding energy. Also 

h(r) = -Pi{r){P l (r) + -Q 2 (r))-^Q 2 {r)[P 1 (r)+Q A (r)/r^] 

r Sr 

where P\ , Q2 , and Q4 are defined in Eq. (j3J) . 

Keeping only the Coulomb and Uehling potentials, one finds 

2a 

p i( r ) = ~oiZ—(2m e )xo{2m e r) 
/ 2a 

Q 2 (r) = aZ ( 1 + — \xx{2m e r) + (2m e r)xo(2m e r)] 
„ , s 2a f°° (z 2 - 1)V2 / 1 

(6qr — (gr) 3 ) cos(gr) + (3(gr) 2 — 6) sin(gr) ^ 



n J Jo I 2 + q 

where Xn{x) is defined in pQ. 

Since vacuum polarization is assumed to be a relatively small correction to the Coulomb 
potential, it will be sufficient to approximate Q2(r) by aZ/r. After some algebra, one can 
reduce the expectation values to single integrals: 



(az) 2 — az + 1 f 



'1 , ^5 ^° + TT, \5 Stt dz 

v f + az)° (1 + az)° J 

When multiplied by —2Bo/(m^ + mjv) this results in a shift of -0.00015 meV for the 2s-state 
and of -0.00001 meV for the 2p-state. For muonic deuterium, the corresponding numbers are 
-0.000176 meV and -0.000030 meV, respectively. 



,aZ , , s o 2a f°° (z 2 - I) 1 / 2 ( 1 

r 07T Vi z \ 2z z 

'2(a^) 2 + l r 1 \ , 

2(1 + az) 4 2(1 + az) 4 / 

When multiplied by \/{m^ + mjv) this results in a shift of 0.00489 meV for the 2s-state and 
of 0.00017 meV for the 2p-state of muonic hydrogen. For muonic deuterium, the corresponding 
numbers are 0.005543 meV and 0.000206 meV, respectively. 

These expectation values also appear when vacuum polarization is included in the Breit 
equation. 
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Finally, 



<f^4(r)> 



{aZfm 2 2a 

37T ./ 1 



6 



oo (z 2 



1)1/2 



1 + 



2z 2 



6 ,2 + az 



2 . 3(az) 2 + 2az- 1 

- In 1 + az) + V / 3 + 

+ az az (1 + az) 



3 + az 



4(1 + az) 



1 - 3az - 2(az) 2 
4(1 + az) 4 



dz 



When multiplied by l/(mu + mjy) this results in a shift of 0.002475 meV for the 2s-state 
and of 0.000238 meV for the 2p-state. For muonic deuterium, the corresponding numbers are 
0.002753 meV and 0.000281 meV, respectively. 

Combining these expectation values according to equations |2 and |21 one finds a contribution 
to the 2p-2s transition of -0.00419 meV (hydrogen) and -0.00479 meV (deuterium). To obtain the 
full relativistic and recoil corrections, one must add the difference between the expectation values 
of the Uehling potential calculated with relativistic and nonrelativistic wave functions, giving a 
total correction of 0.0166 meV for muonic hydrogen. This is in quite good agreement with the 
correction of .0169 meV calculated by Veitia and Pachucki [HHj. The treatment presented here 
has the advantage of avoiding second order perturbation theory. For deuterium, one obtains a 
total correction of 0.0179 meV. 



Appendix 2: Details of Corrections to the Hyperfine Structure of the 2s-state 
of Muonic Hydrogen and Deuterium 

The expectation value of V SlS2 in an ns state with j = 1/2 is 

. „ 2uoa(aZm r ) 3 , ^ r „,„ . , , , 

AE <- - = • (i + "" >[f(f +r >- s ^ + :> - 3/41 

As was shown in 12] , the energy shift of the 2s-state has to be multiplied by: 

1 + &VP + ^vertex + ^Breit + ^FS,rec 

Here ([22]) 



and (0, Eq. (277)) 



2a(aZ) f, , , 13\ , 

Evertex = 3 ' ( Ill(2) - — J = -1.36 ■ 10~ 4 



_ 17(aZ) 2 _ , 

^Breit — 7, — l-lo • ±U 



The vacuum polarization correction has two contributions. One of these is a result of a modifi- 
cation of the magnetic interaction between the muon and the nucleus and is given by (see |16j ^ 



£vpi= ^l r2dr {Wi) I ^ r)GMdq (12) 



(z 2 -!) 1 / 2 f 1+ 1\ dz 



z 2 \ 2z 2 J im 2 [z 2 + (g/2m e ) 2 ] 
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One can do two of the integrals analytically and obtains for the 2s-state (with a = 2m e /(aZm r 
and sinh(0) = q/(2m e ) = K/a) 



K 2 



4a 



3vr 2 J (1 + K 2 ) 2 



F(4>)Gm (aZm r K) dK 



+ 



6 



(1 + K 2 ) (l + K 2 ) 2 



(13) 



where F((p) is known from the Fourier transform of the Uehling potential (given as U2(q) in 
Ref. pQ) and is given by 

1 37T 

F(<t>) = o + (coth 2 (0) - 3) • [1 + <f> ■ coth(^)] = — U 2 (q) 
6 a 

with sinh((/>) = q/2m e . 

The other contribution, as discussed by |38l 13*^] arises from the fact that the lower energy 
hyperfme state, being more tightly bound, has a higher probability of being in a region where 
vacuum polarization is large. This results in an additional energy shift of 

2 / V Uehl (r)f 2s {r)SMip2s(r)d 3 r 



Following Ref. 38 with y = (aZm r /2) ■ r, one has 



exp(-y) 



(l-y)(ln(2y)+ 7 ) + 



aZm r / 

(7 is Euler's constant), and 

^2s(r) = V2s(0)(l - y)exp(-y) 
One finds after a lengthy integration 
16a f 00 dK 



13y - 3 - 2y 2 1 
4 



£VP2 



3tt 2 J l + K 2 



G E (aZm r K)F(<t>) 



1 



17 



+ 



41 



24 



+ 



2 (l + K 2 ) 2 (l + K 2 ) 3 (l + K 2 ) 4 
ln(l + K 2 



l + K 2 



6 



[l + K 2 ) (l + K 2 ) 2 



+ 



tan' 1 (if) 
K 
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+ 



20 



12 



2(1 + K 2 ) (l+K 2 ) 2 (l + K 2 f 



(14) 



Sternheim 39 denotes the two contributions by 5m an d 5e, respectively. An alternative exres- 
sion, obtained by assuming a point nucleus, using Eq.(131) from 1 for the Uehling potential, 
and doing the integrations in a different order, is 



£VP2 



16a 
~3V 



(Z 2 - 1)V2 



1 



az 



+ 



1 

l+ 2^) (l + az) 2 
23 3 



2 l + az 8(1 + az) 2 2(1 + az) 3 



(15) 



+ ln(l + az) ■ 1 



+ 



l + az 2(l + az) : 



dz 
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with a = 2m e / {aZm re d)- Both methods give the same result. 

In the case of ordinary hydrogen, each of these contributes 3a 2 /8 = 1.997 • 10~ 5 . The accuracy 
of the numerical integration was checked by reproducing these results. One can thus expect that 
muonic vacuum polarization will contribute 3a 2 /4 ~ 4- 10 -5 , as in the case of normal hydrogen. 
This amounts to an energy shift of 0.0009 meV in muonic hydrogen and 0.0002 meV in muonic 
deuterium. Contributions due to the weak interaction or hadronic vacuum polarization should 
be even smaller. For muonic hydrogen, one obtains ev / Pi=0-00211 and eyp2=0. 00325 for a 
point nucleus. Including the effect of the proton size (with Ge(q) = Gm(q) as a dipole form 
factor) reduces these numbers to 0.00206 and 0.00321, respectively. For the case of muonic deu- 
terium, the corresponding numbers are ey P i=0.00218 (0.00207) and £^2=0.00337 (0.00326), 
respectively. The contribution to the hyperfine splitting of the 2s-state of hydrogen is then 
0.0470 meV+0. 0733 meV=0. 1203 meV (0.1212 meV if muonic vacuum polarization is included). 
The combined Breit and vertex corrections reduce this value to 0.1207 meV. (0.1226 meV if the 
proton form factors are not taken into account). 

The contribution to the hyperfine structure from the two loop diagrams JH] can be calculated 
by replacing U2{aZm r K) = (a/3ir)F(<p) by Ui(aZm r K) (as given in in equations IT31 and 

1141 The resulting contributions are 1.64 • 10 -5 and 2.46 • 10 -5 (for deuterium 1.69 • 10 -5 and 
2.54 • 10 -5 ), respectively, giving a total shift of 0.0009 meV in muonic hydrogen and 0.0002 meV 
in muonic deuterium. 
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